rl multiplicity and the Witten-Veneziano relation at finite temperature 
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We discuss and propose the minimal generalization of the Witten-Veneziano relation to finite tem- 
peratures, prompted by STAR and PHENIX experimental results on the multiplicity of ri' mesons. 
After explaining why these results show that the zero-temperature Witten-Veneziano relation can- 
not be straightforwardly extended to temperatures T too close to the chiral restoration temperature 
Tch and beyond, we find the quantity which should replace, at T > 0, the Yang-Mills topological 
susceptibility appearing in the T = Q Witten-Veneziano relation, in order to avoid the conflict with 
experiment at T > 0. This is illustrated through concrete T-dependences of pseudoscalar meson 
masses in a chirally well-behaved, Dyson-Schwinger approach, but our results and conclusions are 
of a more general nature and, essentially, model-independent. 
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I. INTRODUCTION 

The ultrarclativistic heavy-ion collider facilities like 
RHIC at BNL and LHC at CERN strive to produce a 
new form of hot QCD matter. The experiments show 
P, Q that it has very intricate properties and presents 
a big challenge especially for theoretical understanding. 
While above the (pseudo) critical temperature Tc ^ 170 
MeV this matter is often called the quark-gluon plasma 
(QGP), it cannot be a perturbatively interacting quark- 
gluon gas (as widely expected before RHIC results 
until significantly higher temperatures T 3> Tc. Instead, 
the interactions and correlations in the hot QCD matter 
are still strong (e.g., see Refs. 0, Q) so that its more 
recent and more precise name is strongly coupled QGP 
(sQGP) One of its peculiarities seems to be that 
strong correlations in the form of quark-antiquark [qq) 
bound states and resonances still exist 1, 13 in thesQGP 
well above Tc- In the old QGP paradigm, even deeply 
bound charmonium (cc) states such as J/'^ and rjc were 
expected to unbind at T « Tc, but lattice QCD simula- 
tions of mesonic correlators now indicate they persist till 
around 2Tc 0,01 or even above Q. Similar indications for 
light-quark mesonic bound states are also accumulating 
from lattice QCD Q and from other methods [1, [13, [HI ■ 
This agrees well with the findings on the lattice (e.g., see 
Ref. [12| for a review) that for realistic explicit chiral 
symmetry breaking (ChSB), i.e., for the physical values 
of the current quark masses, the transition between the 
hadron phase and the phase dominated by quarks and 
gluons, is not an abrupt, singular phase transition but 
a smooth, analytic crossover around the psewdocritical 
temperature Tc- It is thus not too surprising that a clear 
experimental signal of, e.g., deconfinement, is still hard 
to find and identify unambiguously. 

The most compelling signal for production of a new 



form of QCD matter, i.e., sQGP, would be a restora- 
tion - in hot and/or dense matter - of the symmetries 
of the QCD Lagrangian which are broken in the vac- 
uum. One of them is the [SUyi(-^/) flavor] chiral symme- 
try, whose dynamical breaking results in light, (almost- 
)Goldstone pseudoscalar (P) mesons - namely the octet 
P = tt", 71=*=, i4r°, A'", K^,r], as we consider all three light- 
quark flavors, Nf = 3. The second one is the Ua(1) 
symmetry. Its breaking by the non-Abelian axial Adler- 
Bell-Jackiw anomaly ('gluon anomaly' for short) makes 
the remaining pseudoscalar meson of the light-quark sec- 
tor, the r/', much heavier, preventing its appearance as 
the ninth (almost-)Goldstone boson of dynamical chiral 
symmetry breaking (DChSB) in QCD. 

The first experimental signature of a partial restora- 
tion of the Uyi(l) symmetry seems to have been found 
in the ^/snn = 200 GeV central Au-|-Au reactions at 
RHIC. Namely, Csorgo et al. 1^ 14| analyzed combined 
data of PHENIX [H] and STAR [l^ collaborations very 
robustly, through six popular models for hadron multi- 
plicities, and found that at 99.9% confidence level, the 
T]' mass, which in the vacuum is M^/ = 957.8 MeV, is 
reduced by at least 200 MeV inside the fireball. It is the 
sign of the disappearing contribution of the gluon axial 
anomaly to the 77' mass, which would drop to a value 
readily understood together with the (flavor-symmetry- 
broken) octet of qq' (g, q' = u, d, s) pseudoscalar mesons. 
This is the issue of the "return of the prodigal Goldstone 
boson" predicted [l3| as a signal of the Ua(1) symmetry 
restoration. 

Another related but less obvious issue to which we 
want to draw attention, concerns the status, at T > 0, of 
the famous Witten-Veneziano relation (WVR) [isl [H 
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between the 77', rj and ii'-meson masses M,ji^ri,K, pion 
decay constant , and Yang-Mills (YM) topological sus- 
ceptibility xym • WVR was obtained in the limit of large 
number of colors Nc [l^, [l^ . It is well satisfied at T = 
for Xym obtained by lattice calculations (e.g., (20l - [23| ). 
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Nevertheless, the T-dependence of xym is such that 
the straightforward extension of Eq. (P) to T > [23 |. 
i.e., replacement of all quantities^ therein by their re- 
spective T-dependent versions Mr,'{T), M^{T), Mk{T), 
f-n{T) and xym{T), leads to a conflict with experiment 
|13l 111] . Since this extension of Eq. dJ) to T > was 
studied in Ref. [l^] before the pertinent experimental 
analysis [l^, [3 > O'^^ of the purposes of this paper is to 
revisit the implications of the results of Ref. [2^ for 
WVR at T > 0, and demonstrate explicitly that they are 
practically model-independent. The other, more impor- 
tant purpose is to propose a mechanism which can enable 
WVR to agree with experiment at T > 0. 



II. THE RELATIONS CONNECTING TWO 
THEORIES, QCD AND YM 

Both issues pointed out before Eq. ^ and around it, are 
best understood in a model-independent way if one starts 
from the chiral limit of vanishing current quark masses 
{vriq = 0) for all three light flavors, q = u,d, s. Then not 
only pious and kaons are massless, but is also rj, which is 
then (since the situation is also SU(3)-flavor-symmetric) 
a purely SU(3)-octet state, rj = rjg. In contrast, r/ is then 
purely singlet, rj' = r]o; since the divergence of the singlet 
axial quark current 97^75 ^A'^q is nonvanishing even for 
niq = due to the gluon anomaly, the r/' mass squared 
receives the anomalous contribution AM^, A"*//^, in 
the notation of Ref. [iTj ) which is nonvanishing even in 
the chiral limit: 

|-=AM,^„=AA4.^+0(i-). (2) 

However, A'' and /,,' are known accurately^ only in the 
large limit. There, in the leading order in l/iVc, A^ 
is given by the YM (i.e., "pure glue" ) topological suscep- 
tibility Xym times 2Nf = 6 llSl. [l9ll . and the "77' decay 
constant" is the same as /^r Thus, keeping only 
the leading order in 1/Nc, the last equality is WVR in 
the chiral limit. 

The consequences of Eq. ([2]) remain qualitatively the 
same realistically away from the chiral limit. This will 
soon become clear on the basis of, e.g., Eq. ([S]) below. 
Namely, due to DChSB in QCD, for relatively light cur- 
rent quark masses niq {q = u,d,s), the qq' bound-state 
pseudoscalar meson masses (including the nonanomalous 



parts of the 77' and 77 masses) behave as 

M|g, = const {niq + niqi) , (q, q' = it, d, s). (3) 

The pseudoscalar mesons (including ?/) thus obtain rel- 
atively light nonanomalous contributions Mqqi to their 
masses Mp, allowing them to reach the empirical val- 
ues. That is, instead of the eight strictly massless Gold- 
stone bosons, 7rO,7r±,i^o,i^o,if± and 77 are relatively 
light almost-Goldstones. Among them, in the limit of 
isospin symmetry {niu ~ rrid), only rj now receives also 
the gluon-anomaly contribution since the explicit SU(3) 
flavor breaking between the nonstrange (NS) u, d-quarks 
and s-quarks causes the mixing between the isoscalars 77 
and 77'. For rriq ^ 0, Eq. ([2]) is replaced by the usual 
WVR ^ containing also the nonanomalous contribu- 
tions to meson masses. Nevertheless, these contributions 
largely cancel due to the approximate SU(3) flavor sym- 
metry and to DChSB [i.e., Eq. ©]. 

This can be seen assuming the usual SU(3) qq con- 
tent of the pseudoscalar meson nonet with well-defined 
isospin'^ quantum numbers, in particular the isoscalar 
(/ = 0) octet and singlet etas, 773 = {uii + dd — 2ss)/-\/6, 
770 = {uu + dd + ss) I >/3, whose mixing yields the phys- 
ical particles 77 and 77'. Since the nonanomalous parts 
of the 770 and 778 masses squared, Af^^ and M^^, are re- 
spectively M^f^ - jA^i + and Affg sa |M|. ~ ^M^ 
(see, e.g., Ref. HSf), and since M^^ + M^^ = Af^ ^ ]^2^ 
the nonanomalous parts of the 77 and rj' masses arc can- 
celed by 2A/|. in WVR (P). Another way of seeing this 
is expressing the nonanomalous parts of A/,^ + Af^, = 
A/2^ + by Eq. Q. Thus again Af^ + M^, - 2M]^ « 
AA/^j^, showing again that already WVR's chiral- limit- 
nonvanishing part ([2]) reveals the essence of the influence 
of the gluon anomaly on the masses in the 77'-?; complex. 
This is important also for the presently pertinent finite- 
T context because thanks to this, below it will be shown 
model-independently that WVR ([1]) containing the YM 
topological susceptibility xym implies T-depcndence of 
7/ mass in conflict with the recent experimental results 

MM- 

Namely, the gluon-anomaly contribution ([2|) is estab- 
lished at r = but it is not expected to persist at 
high temperatures. Ultimately, 77' should also become 
a massless Goldstone boson at sufficiently high T, where 
Xym{T) 0. However, according to WVR, AAf^/(r) 
falls only for T where /^(T)^ does not fall faster than 
^Xym(T), as stressed in Ref. p^ . 

The WVR's chiral-limit version ^ manifestly points 
out the ratio Xym{T)/ fTr{T)^ as crucial for the anomalous 



^ Throughout this paper, all quantities are for dofinitoness as- 
sumed at T = unless their T-dependence is specifically in- 
dicated in formulas or in the text. 

^ Also note that a unique decay constant" f^i is, strictly speak- 
ing, not a well-defined quantity, as two rj' decay constants are ac- 
tually needed: the singlet one, , and the octet one, ; e.g., 
see an extensive review [25II or the short Appendix of Ref. [2^ . 



^ Vanishing of the anomaly at sufficiently high T opens the pos- 
sibility of studying the interesting scenario of maximal isospin 
violation at high T [28l . |29|| , but as the effects of the small differ- 
ence between rriu and are not important for the present con- 
siderations, we stick to the isospin limit throughout the present 
paper. 
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T]' mass, but the above discussion shows that this remains 
esscntiahy the same away from the chiral hmit. 

In the present context, it is important for practical 
calculations to go realistically away from the chiral limit, 
in which the chiral restoration is a sharp phase transi- 
tion at its critical temperature Tch where the chiral-limit 
pion decay constant vanishes very steeply, i.e., as steeply 
as the chiral quark condensate. In contrast, for realis- 
tic explicit ChSB, i.e., m„ and of several MeV, this 
transition is a smooth crossover (e.g., see Ref. [l^). For 
the pion decay constant, this implies that fTr{T) still falls 
relatively steeply around pseudocritical temperature Tch, 
but less so than in the chiral case, and even remains fi- 
nite, enabling the usage of WVR ([T]) for the temperatures 
across the chiral and Ua(1) symmetry restorations. 

WVR is very remarkable because it connects two dif- 
ferent theories: QCD with quarks and its pure-gauge, 
YM counterpart. The latter, however, has much higher 
characteristic temperatures than QCD with quarks: the 
"melting temperature" Tym where xym(2^) starts to de- 
crease appreciably was found on lattice to be, for exam- 
ple, Tym~ 260 MeV [H [H or even higher, Tym « 300 
MeV [331 • In contrast, the pseudocritical temperatures 
for the chiral and deconfincment transitions in the full 
QCD are lower than Tym by some 100 McV or more (e.g., 
see Ref. [l^l ) due to the presence of the quark degrees of 
freedom. 

This difference in characteristic temperatures, in con- 
junction with xym{T) in WVRs ^ and ^ would im- 
ply that the (partial) restoration of the 11^1(1) symme- 
try (understood as the disappearance of the anomalous 
770/7/ mass) should happen well after the restoration of 
the chiral symmetry. But, this contradicts the RHIC ex- 
perimental observations of the reduced if mass [l^ [T^ | 
i/WVRs ([!]), ^ hold unchanged also close to the QCD 
chiral restoration temperature T ny,, around which fj^iT) 
decreases still relatively steeply^ [2J| for realistic explicit 
ChSB, thus leading to the increase of Qxym{T) / fyriT)^ 
and consequently also of Af^'. 

There is still more to the relatively high resistance of 
Xym(^) to temperature: not only does it start falling 
at rather high Tym, but xym(T) found on the lattice 
is falling with T relatively slowly. In some of the ap- 
plications in the past (e.g., see Refs. [U, H^), it was 
customary to simply rescale a temperature characteriz- 
ing the pure-gauge, YM sector to a value characterizing 
QCD with quarks. (For example, Refs. [sil . Issj rescaled 
Tym = 260 McV found by Ref. [3l| to 150 MeV). How- 
ever, even if we rescale the critical temperature for melt- 
ing of the topological susceptibility XYuiT) from Tym 
down to Tch, the value of 6xYM(7')//7r(T)^ still increases 
a lot [23 | for the pertinent temperature interval starting 



already below Tch- This happens because yym(T) falls 
with T more slowly than /7r(T)^. (It was found [2J] that 
the rescaling of Tym would have to be totally unrealistic, 
to less than 70% of Tch, in order to achieve sufficiently 
fast drop of the anomalous contribution that would allow 
the observed enhancement in the 77' multiplicity.) 

These WVR-induced enhancements of the 77' mass for 
T ^ Tch were first noticed in Ref. [2^. This reference 
used a concrete dynamical model (with an effective, rank- 
2 separable interaction, convenient for computations at 
T > 0) HI] of low energy, nonperturbative QCD to ob- 
tain mesons a s gg ' bound states in Dyson-Schwinger (DS) 
approach [STHSOj , which is a bound-state approach with 
the correct chiral behavior ([3]) of QCD. Nevertheless, this 



concrete dynamical DS model was used in Ref. [2J| to 
get concrete values for only the nonanomalous parts of 
the meson masses, but was essentially not used to get 
model predictions for the mass contributions from the 
gluon anomaly, in particular xym{T). On the contrary, 
the anomalous mass contribution was included, in the 
spirit of l/Nc expansion, through WVRJT]). Thus, the 
T-evolution of the 77'-77 complex in Ref. [2J] was not dom- 
inated by dynamical model details, but by WVR, i.e., 
the ratio 6xym(T)//^(T)2. Admittedly, /^(T) was also 
calculated within this model, causing some quantitative 
model dependence of the anomalous mass in WVR, but 
this cannot change the qualitative observations of Ref. 
[24| on the 77' mass enhancement. Namely, our model 
/7r(T), depicted as the dash-dotted curve in Fig. [U obvi- 
ously has the right crossover features [l2| . It also agrees 
qualitatively with fviTys calculated in other realistic dy- 
namical models [l0|,|38|. Various modifications were tried 
in Ref. 2J] but could not reduce much the 77' mass en- 



hancement caused by this ratio, let alone bring about the 
significant 7/ mass reduction found in the RHIC experi- 
ments [Hill. 

One must therefore conclude that either WVR breaks 
down as soon as T approaches Tch, or that the T- 
dependence of its anomalous contribution is different 
from the pure-gauge xym{T). We will show that the 
latter alternative is possible, since WVR can be recon- 
ciled with experiment thanks to the existence of another 
relation which, similarly to WVR, connects the YM the- 
ory with full QCD. Namely, using large-A''c arguments, 
Leutwyler and Smilga derived [l^, at T = 0, 



Xym 



X 



1 I V 

A rn {qq)o 



(4) 



the relation (in our notation) between the YM topologi- 
cal susceptibility xym, and the full-QCD topological sus- 
ceptibility Xj the chiral-limit quark condensate {qq)o, and 
TO, the harmonic average of Nj current quark masses m^. 
That is, TO is Nf times the reduced mass. In the present 
case of Nf = 3, q = u, d, s, so that 



* Relative to decay constants of mesons containing a strange quark; 
e.g., compare fss{T) of the unphysical ss pseudoscalar with 
U{T) in Fig.[T] 



^ = y — . 

m ^ jjiq 

q—u.d,s 



(5) 
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Equation (|4]) is a remarkable relation between the two 
pertinent theories. For example, in the limit of all very 
heavy quarks (m^ — > oo, q = w, d, s), it correctly leads to 
the result that xym is equal to the value of the topologi- 
cal susceptibility in quenched QCD, xym = x(TOi} = oo). 
This holds because x is by definition the vacuum expec- 
tation value of a gluonic operator, so that the absence of 
quark loops would leave only the pure-gauge, YM con- 
tribution. However, the Leutwyler-Smilga relation ([4]) 
also holds in the opposite (and presently pertinent) limit 
of light quarks. This limit still presents a problem for 
getting the full-QCD topological susceptibility x on the 
lattice [i^l, but we can use the light-quark-sector result 

MM 
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N 
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(6) 



where Cm stands for corrections of higher orders in small 
niq, and thus of small magnitude. The leading term is 
positive (as {qq)o < 0), but Cm is negative, since Eq. dH) 
shows that x < min(— to {qq)o/Nf,XYM)- 

Although small, C,„ should not be neglected, since 
Cm = would imply, through Eq. that xym = oo. 
Instead, its value (at T = 0) is fixed by Eq. (|4]): 



Cm ^?7i(0) 



t (gg)c 
Nf 



1 — Xym 
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(7) 



All this starting from Eq. ([4]) has so far been at T = 0. 
If the left- and right-hand side of Eq. ([4]) are extended to 
T > 0, it is obvious that the equality cannot hold at ar- 
bitrary temperature T > 0. The relation (|4| must break 
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FIG. 1: The relative-temperature dependences, on T/Tch, of 
X^^'^, {qq)y^ , frr and fss, i.e., the T/Tch-dependences of the 
quantities entering in the anomalous contributions to various 
masses in the i]'-ri complex - see Eq. (|10|l and formulas below 
it. The solid curve depicts x''^^ for 5 = in Eq. (|9}, and 
the short-dashed curve is X^''* for S = 1. At T = 0, the both 
X 's are equal to xym = (0.1757 GeV)*, the weighted average 
[24| of various lattice results for xym . The dotted (red) curve 
depicts ~{qq)]J^ , the dash-dotted (blue) curve is and the 
long-dashed (blue) curve is fss- (Colors online.) 



down somewhere close to the (pseudo)critical tempera- 
tures of full QCD (~ Tch) since the pure-gauge quantity 
Xym is much more temperature-resistant than the right- 
hand side, abbreviated as x- The quantity x, which may 
be called the effective susceptibility, consists of the full- 
QCD quantities x and {qq)o, the quantities of full QCD 
with quarks, characterized by Tch, just as f-n{T). As 
T Tchj the chiral quark condensate {qq)Q(T) drops 
faster than the other DChSB parameter in the present 
problem, namely firiT) for realistically small explicit 
ChSB. (See Fig. [T]for the results of the dynamical model 
adopted here from Ref. [H, and, e.g., Refs. [H, Ull 
for analogous results of different DS models). Thus, the 
troublesome mismatch in T-dependences of fniT) and 
the pure- gauge quantity xym(T'), which causes the con- 
flict of the temperature-extended WVR with experiment 
at T ^ Tch, is expected to disappear if xym{T) is re- 
placed by x(2^)i the temperature-extended effective sus- 
ceptibility. The successful zero-temperature WVR ([T]) is, 
however, retained, since xym = x at T = 0. 

Extending Eq. ([6]) to T > is something of a guess- 
work as there is no guidance from the lattice for x{T) [un- 
like Xym (7^)]- Admittedly, the leading term is straight- 
forward as it is plausible that its T-dependence will sim- 
ply be that of {qq)o{T). Nevertheless, for the correction 
term Cm such a plausible assumption about the form of 
T-dependence cannot be made and Eq. ([7]), which re- 
lates YM and QCD quantities, only gives its value at 
T = 0. We will therefore explore the T-dependence of 
the anomalous masses using the following Ansatz for the 
T > generalization of Eq. ([6]): 



x(t) = - 



m{qq)Q{T) 
Nf 



CmiO) 



{qq)oiT) 



(gg)o(T = 0) 



(8) 



where the correction-term T-dependence is parametrized 
through the power S of the presently fastest- vanishing (as 
T -> Tch) chiral order parameter {qq)o{T). 

The T > extension ^ of the light-quark x, Eq. ([6|, 
leads to the T > extension of x- 



X(T) = 



m{qq)o{T) 



TV. 



/ 



m{qq)o{T) 

NfCmiO) 



{qq)o{T = 0) 



-I <5^ 



(qqMT) 



(9) 

We now use x(^) iu WVR instead of xym(T) used by 
Ref. [l^]. This gives us the temperature dependences of 
the masses in the r]-r]' complex, such as those in Fig. [2] 
iUustrating the cases i5 = and S ~ 1. 

It is clear that x(T) Q blows up as T — > Tch if the cor- 
rection term there vanishes faster than {qq)o{T) squared. 
Thus, varying S between and 2 covers the cases from the 
T-independent correction term, to (already experimen- 
tally excluded) enhanced anomalous masses for d notice- 
ably above 1, to even sharper mass blow-ups for 5 — > 2 
when T — > Tch- On the other hand, it does not seem 
natural that the correction term vanishes faster than the 
fastest- vanishing order parameter (gg)o(T). Indeed, al- 
ready for the same rate of vanishing of the both terms 
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{6 = 1), one can notice in Fig. [2] the start of the precur- 
sors of the blow-up of various masses in the rj'-r] complex 
as T — > Tch although these small mass bumps are still 
experimentally acceptable. Thus, in Fig. [2]we depict the 
6=1 case, and 6 = (T-independent correction term) as 
the other acceptable extreme. Since they turn out to be 
not only qualitatively, but also quantitatively so similar 
that the present era experiments cannot discriminate be- 
tween them, there is no need to present any 'in-between 
results', for < (5 < 1. 

Next we turn to completing the explanation how the 
above-mentioned results in Fig. [5] were obtained. 
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FIG. 2: The relative-temperature dependence, on T/Tch, of 
the pseudoscalar meson masses for two x{T), namely Eq. (|9]) 
with 5 = (upper panel) and with S = 1 (lower panel). The 
meaning of all symbols is the same on the both panels: the 
masses of r;' and rj are, respectively, the upper and lower solid 
curve, those of the pion and nonanomalous ss pseudoscalar 
are, respectively, the lower and upper dash-dotted curve, MrjQ 
and Mr,^ are, respectively, the short-dashed (red) and long- 
dashed (red) curve, Mn^g is the medium-dashed (blue), and 
M,,g is the dotted (blue) curve. (Colors online.) The straight 
line 2-kT is twice the lowest Matsubara frequency. 



Using x{T) in WVR instead of Xym{T) used by Ref. 
[23 |. does not change anything at T = 0, where x(T) = 
Xym(O), which remains an excellent approximation even 
well beyond T = 0. Nevertheless, this changes drasti- 
cally as T approaches Tch- For T ~ Tch, the behavior 
of x{T) is dominated by the T-dependence of the chiral 
condensate, tying the restoration of the \]a{1) symmetry 
to the chiral symmetry restoration. 

As for the nonanomalous contributions to the me- 
son masses, we use the same DS model (and param- 
eter values) as in Ref. since it includes both 
DChSB and correct QCD chiral behavior as well as re- 
alistic explicit ChSB. That is, all nonanomalous results 
(A/tt, /it, A^if , /if, as well as Mss and fss, the mass and 
decay constant of the unphysical ss pseudoscalar meson) 
in the present paper are, for all T, taken over from Ref. 
2J] . We used this same model also for computing the chi- 
ral quark condensate {qq)o, including its T-dcpendence 
displayed in Fig. [T] 

This defines completely how the results displayed in 
Fig . [2]were generated. For details, see Ref. [131 {^ind Ref. 
}30| for and M^^ ). Here we list only the formulas 
which, in conjunction with Fig. [U enable the reader to 
understand easily the T-dependences of the masses in 
Fig. [5J The theoretical 77' and rj mass eigenvalues are 



A4(r) 

M2(T) 



I [M^.s (T) + Ml (T) + A,,, (T)] , (10) 



m + MliT) 



'(T)], (11) 



where A„ 



Vns 



AP P 



R = ^ ^ Y = !jL- 



aP 

VNS 



Mi 



M^ + 2(3, 



nP 



AI% + (3X' 



-i2 + Xfp,M, 



AL 



88 + -(1-X)^/3, 



AL 



88 



_ M^- -I- - AP 



AI, 



00 



\ AI% -f ^ All 
(|lip . the T-depcndencc is 



In all expressions after Eq, 
understood. 

In both cases considered for the topological suscepti- 
bility ([5]) [(5 = 0, i.e., the constant correction term, and 
(5 = 1, i.e., the strong T-dependence cx {qq)o{T) of both 
the leading and correction terms in x(^)]: the results are 
consistent with the experimental findings on the decrease 
of the rj' mass of Csorgo et al. fl^]. 



III. SUMMARY, DISCUSSION AND 
CONCLUSIONS 



In the light of the recent experimental results on the 
rj' multiplicity in heavy- ion collisions [l3l [l^ , we revis- 
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ited the earlier theoretical work [23| concerning the ther- 
mal behavior of the ry'-ry complex following from WVR 
straightforwardly extended to T > 0. We have con- 
firmed the results of Ref. [l^l on WVR where the ratio 
Xym{T)/ fTr{T)^ dominates the T-dependence, and clari- 
fied that these results are practically model- independent. 
It is important to note the difference between our ap- 
proach and those that attempt to give model predic- 
tions for topological susceptibility, such as Refs. [42|.|43|. 
By contrast, in Refs. [2J] and here, as well as earlier 
works m, ia, S, m at T = 0, a DS dynamical model 
is used (as far as masses are concerned) to obtain only 
the nonanomalous part of the light pseudoscalar meson 
masses (where the model dependence is however dom- 
inated by their almost- Goldstone character), while the 
anomalous part of the masses in the ry'-ry complex is, 
through WVR, dictated by 6xym//^- In this ratio, /7r(T') 
is admittedly model-dependent in quantitative sense, but 
other realistic models yield qualitatively similar crossover 
behaviors [i^ of f-n{T) for 0, as exemplified by our 

pg. m and Fig. 2 in Ref. [Hi, and by Fig. 6 in Ref. 
p^ . Such /^(T) behaviors are also in agreement with 
the T-dependence expected of the DChSB order param- 
eter on general grounds: a pronounced fall-off around 
Tc\i ~ but exhibiting, in agreement with lattice [l^l, a 
smooth crossover pattern for nonvanishing explicit ChSB, 
a crossover which gets slower with growing niq [e.g., com- 
pare /^(T) with fss{T) in Fig. [T]. In contrast to the QCD 
topological susceptibility x, the YM topological suscep- 
tibility and even its T-dependence Xym{T), including its 
"melting" temperature Tym, can be extracted [l^l rea- 
sonably reliably from the lattice [H,|3l|. Thus, it was not 
modeled in Ref. [IJ] . Hence our assertion that the results 
of Ref. [iSl unavoidably imply that the straightforward 
extension of WVR to T > is falsified by experiment 
[13I especially if one recalls that even the sizeable 
T-rescaling [sl, [l^l Tym — > Tch was among the attempts 
to control the rj' mass enhancement [2^ . 

Nevertheless, we have also shown that there is a plau- 
sible way to avoid these problems of the straightforward, 
naive extension of WVR to T > 0, and this is the main re- 
sult of the present paper. Thanks to the existence of an- 
other relation, Eq. ([4]), connecting the YM quantity xym 
with QCD quantities x ^-nd {qq)o, it is possible to define 



a quantity, x, which can meaningfully replace xym(^) in 
finitc-T WVR. It remains practically equal to xym up to 
some 70% of Tch, but beyond this, it changes following 
the T-dependence of {qq)o{T). In this way, the success- 
ful zero-temperature WVR is retained, but the (partial) 
restoration of Uyi(l) symmetry [understood as the disap- 
pearing contribution of the giuon anomaly to the rj' (770) 
mass] is naturally tied to the restoration of the SUyi(3) 
flavor chiral symmetry and to its characteristic temper- 
ature Tch, instead of Tym- 

It is very pleasing that this fits in nicely with the recent 
ab initio theoretical analysis using functional methods 
j47| . which finds that the anomalous breaking of Uyi(l) 
symmetry is related to DChSB (and confinement) in a 
self-consistent manner, so that one cannot have one of 
these phenomena without the other. 

Of course, the most important thing is that this ver- 
sion of the finite-T WVR, obtainedby xym(T) ^ x{T), 
is consistent with experiment [l^ [l3l for all reasonable 
strengths of T-dependence [0 < 5 < 1 in Eq. ©]. 
Namely, the both tablets in Fig. [2] show, first, that rj' 
mass close to Tch suffers the drop of more than 200 MeV 
with respect to its vacuum value. This satisfies the min- 
imal experimental requirement abundantly. Second, Fig. 
[5] shows an even larger drop of the 770 mass, to some 400 
MeV, close to the "best" value of the in-medium 77' mass 
(340 MeV, albeit with large errors) obtained by Csorgo 
et al. [H, . This should be noted because the 770 mass 
inside the fireball is possibly even more relevant. Namely, 
although it is, strictly speaking, not a physical meson, 770 
is the state with the qq content closest to the qq con- 
tent of the physical 77' in the vacuum. Thus, among the 
isoscalar qq states inside the fireball, 770 has the largest 
projection on, and thus the largest amplitude to evolve, 
by fireball dissipation, into an 77' in the vacuum. 
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